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Abstract. This constribution collects new recent results on preheating after inflation. 
We discuss tachyonic preheating in the SUSY motivated hybrid inflation; development 
of equilibrium after preheating; theory of fermionic preheating and the problem of 
gravitino overproduction from preheating. 

1 Introduction 

The best-fit universe is the flat FRW universe with scale-free scalar metric fluctu- 
ations. The model is in agreement with the observations of the cosmic microwave 
background radiation (CMB) anisotropy, the large scale structure (LSS) of the 
universe and the tests of the global geometry such as the high redshift supernovae 
(SNc). The theoretical scenario of the expanding universe includes a very early 
stage of inflation, subsequent stage of particle creation during preheating and the 
setting of the thermal equilibrium. The resulting hot expanding universe passes 
through a sequence of crucial stages. As the universe expands, it opens up for 
the observational tests, beginning "upward" from the big bang nucleosynthesis 
(BBN), see the left panel of Figure |. 

Although inflation is an essential part of the cosmological paradigm, from 
the particle physics perspective it is not easy to construct a satisfactory model 
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Fig. 1. Sketch of an expanding universe where the wavelengths X(t) are proportional 
to a scale factor a(t) and "reciprocal" universe where the momenta k(t) ~ 1/A(£) are 
inversly proportional to a(t). 
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of the very early universe including inflation. The right panel of Figure |l] shows 
a theoretical "reciprocal" universe where the momenta k ~ 1/A are increasing 
backwards in time. The reciprocal universe at the earliest times at high momenta 
opens up for theoretical possibilities of the early universe physics, including re- 
cent developments in brane cosmology, string motivated cosmology, supergravity 
in cosmology etc. 

Fundamental M-theory should encompass both supergravity and string the- 
ory. At present the low-energy phenomenology is described by the N = 1 d = A 
supergravity. Some preferable choices of the Kahler potentials, superpotentials 
and Yang-Mills couplings hopefully will be selected at the level of the fundamen- 
tal theory. Until the fundamental theory of all interactions is well understood, 
one may try to address the issues of the eary universe cosmology in the context of 
the most general phenomenological N — 1 supergravity- Yang-Mills-matter the- 
ory. This, in fact, was the case during the last almost 20 years. A rather lengthy 
N = 1 phenomenological supergravity Lagrangian begins with the terms 

, 1 9 - • - ^ / WW* 

e- x L = — MlR - du&d^i + e K V l WV t W - 3 T - 

2 \ Mp 

- Xj W - t V>Xi - (f^V^WxiXi + h-c) 

- \%R" + (le^W^Rj^R + Vvl^Vx, + i> R ^e K/2 V l W + h.c\ 

+ - (1) 

A particular choice of the form of the Lagrangian is motivated and notations are 
given in [^J. In Eq. ([!]) K is the Kahler potential, <P % is the complex conjugate of 

According to the inflationary scenario, the Universe initially expands quasi- 
exponentially in a vacuum-like state without entropy or particles. At the stage 
of inflation, all energy is contained in a classical slowly moving fields <P in the 
inflaton sector. The last term of the 1st line of (|l|) is the scalar potential V(4>i). 
The equations of motion based on the first line should describe inflation, which 
is a challenging problem by itself. The Lagrangian ([!]) contains also other fields 
which give subdominant contributions to gravity. The Friedmann equation for 
the scale factor a(t) and the Klein-Gordon equation for @(t) determine the evolu- 
tion of the background fields. In the chaotic inflation models, soon after the end 
of inflation, an almost homogeneous inflaton field <I>(t) coherently oscillates with 
a very large amplitude of the order of the Planck mass around the minimum of its 
potential. This scalar field can be considered as a coherent superposition of infla- 
tons with zero momenta. The amplitude of oscillations gradually decreases not 
only because of the expansion of the universe, but also because energy is trans- 
ferred to particles created by the oscillating field. At this stage we shall recall 
the rest of the fundamental Lagrangian which includes all the fields interacting 
with inflaton. These interactions lead to the creation of many ultra-relativistic 
particles from the inflaton. Gradually, the inflaton field decays and transfers all 
of its energy to the created particles. In this scenario all the matter constituting 
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the universe is created from this process of reheating. If the creation of parti- 
cles is sufficiently slow, the particles would simultaneously interact with each 
other and come to a state of thermal equilibrium at the reheating temperature 
Tr. This gradual reheating can be treated with the perturbative theory of par- 
ticle creation and thermalization. However, typically particle production from 
coherently oscillating inflatons occurs not in the perturbative regime but in the 
non-perturbative regime of preheating H . 

In this contribution I will discuss several problems of preheating after infla- 
tion, some of them related to the supergravity. This discussion is based on new 
results derived in recent papers ||,||,|]|6]|7],§] . 

2 Tachyonic Preheating 

Another popular class of inflationary models - hybrid inflation - involve multiple 
scalar fields <Pi in the inflaton sector. Previous studies of preheating in hybrid 
models were concentrated on particle creation by parametric resonance that may 
occur when homogeneous background fields oscillate around the minimum of the 
potential. Such parametric resonance may or may not be strong depending on 
the coupling parameters. However, we recently found j?]] that there is strong 
preheating in hybrid inflation, but its character is quite different from preheat- 
ing based on parametric resonance. It turns out that typically there is tachyonic 
instability that appears in a broad class of hybrid inflation models. The backre- 
action of rapidly generated fluctuations does not allow homogeneous background 
oscillations to occur because all energy of the oscillating field is transferred to 
the energy of long-wavelength scalar field fluctuations within a single oscillation. 
However, this does not preclude the subsequent decay of the Higgs and infla- 
ton inhomogeneities into other particles, and thus reheating without parametric 
resonance. 

Consider the simple potential for the two-field hybrid inflation is 

V{4>,o) = \{o 2 -v 2 ) 2 + ^4> 2 o 2 , (2) 

where we used notations <Pi — <fi, <&2 — Inflation in this model occurs while 
the homogeneous field slow rolls from large 4> towards the bifurcation point at 
(f> = (due to the slight lift of the potential in direction). Once </)(t) crosses 
the bifurcation point, the curvature of the a field, m 2 = d 2 V/da 2 , becomes 
negative. This negative curvature results in exponential growth of a fluctuations. 
Inflation then ends abruptly in a "waterfall" manner. 

One reason to be interested in hybrid inflation is that it can be implemented 
in supersymmetric theories. In particular, for illustration we will consider pre- 
heating in the supersymmetric F-term inflation as an example of a hybrid model. 

The simplest F-term hybrid inflation model (without undesirable domaine 
walls) is based on a superpotential with three left-chiral superfields <Pi = (<Pi, <£ 2 , ^3) 
in the Lagrangian ([!]) 

W = ^#1 (4<£ 2 <2>3 - v 2 ) . (3) 
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In this case, the spontaneous breaking of the local (global) {7(1) symmetry be- 
tween the and $3 fields will lead to gauge (global) string formation. 

In global SUSY, using the same notation for superfields and their complex 
scalar components, this superpotential contributes 

V F = ^\4$ 2 $3 - v 2 \ 2 + 4\\$i\ 2 (\$ 3 \ 2 + \<P 2 \ 2 ) ■ (4) 

to the scalar potential. In general, <P 3 and <£ 3 could be (oppositely) charged under 
a local f/(l) symmetry, in which case we should include a D-term, Vfj, which we 
negelect here. 

In this model, inflation occurs when chaotic initial conditions lead to (|#i|) ^> 
v. When this happens, the fields ^2 and $3 acquire large effective masses and 
roll to their local minimum at (#2) = (^3) = 0. In this limit, the potential (0) 
becomes V « ^j-, which gives rise to a non- vanishing effective cosmological 
constant. However, this is a false vacuum state; the true vacuum corresponds to 
(^2^3) = x and (#1) = 0. The slow-roll potential drives the evolution of the 
inflaton towards its true VEV. When its magnitude reaches the value (|# c |) = § 
spontaneous symmetry breaking occurs. 

For further discussion of symmetry breaking in this model, let us rewrite (|J) 
in terms of polar fields: $3 = \$z\e , <£ 2 = 1^2 |e* e . The potential becomes 

V F = - (16|<£ 2 | 2 |<£ 3 | 2 - 8v 2 |<Z> 2 ||<P 3 | cos(6 + 8) + v 4 ) + 4A|<£ 2 | 2 (|^ 2 | 2 + |^ 3 | 2 ) ■ 

At the stage of symmetry breaking, when (|#2^3|) begins to move away from 
zero, the absolute phase Arg(^2^3) = 8 + 8 acquires a mass and is forced to 
zero. Note, however, that the potential is independent of the relative phase, 
9 — 8, reflecting the {7(1) symmetry. Thus, in a quasi-homogeneous patch, the 
{7(1) symmetry allows us to choose the relative phase of the ^2, ^3 fields to be 
zero without any loss of generality. This choice, combined with the vanishing of 
the absolute phase, is equivalent to choosing the two complex ^3 fields to 
be real. In order to leave canonical kinetic terms, we define cr± = \^ \ ± \$2\- 
Furthermore, as inflation has left the inflaton homogeneous across all the patches, 
we may choose it to be real: <f) = \/2|^i|- 

In terms of these three real fields, the potential now becomes 

V F = ^(a 2 + -a 2 _ - v 2 f + \cf> 2 (a 2 + + a 2 _) . (6) 

In the symmetric phase, when a± = 0, the a fields have an effective mass- 
squared: m\{<j)) = X(2(j) 2 + v 2 ). We can now see that spontaneous symmetry 
breaking occurs in this model exactly as in the two field model (||). For 4> < 
(f> c = -j=, the (T+ field has a tachyonic mass that triggers symmetry breaking 
and the end of inflation. On the other hand, the cr_ field has always a large and 
positive effective mass-squared, pinning it to zero. Thus, during inflation and at 
the initial stages of symmetry breaking, this model behaves just like the standard 
two field hybrid model discussed in the last section. We have only to apply the 
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Fig. 2. Evolution of the background fields <J+(t) and <p(t) after symmetry breaking and 
the log of V+{k) for the mode with the momentum k — 0.2\/Aw. 

constraint g 2 — 2A and identify the Higgs field with er + . The equations for the 
homogeneous background components <p(t) and a± (i) admits simple solution 

<f>(t) + -j=a+(t) = fa , <r_(t)=0. (7) 

To study preheating in the F-term inflation, we have to analyse evolution 
of the vacuum fluctuations. Consider vacuum fluctuations in the inflaton sector 
#i of the theory (||). Usual description in terms of a homogeneous background 
plus small fluctuations gives us equations for fluctuations around the background 
solution (0). We define the variances of fields: 

<k ± -{<T ± )| 2 > ren = / (0 [\6a k± (t)\ 2 -\6a k± (0)\ 2 ]^ J fp ± (k,t) (8) 

and similar for cj> field. Here T'± y( p(k,t) are the spectra of the fluctuations. 

Numerical solutions of the equations for the background fields and for the 
time evolution of the fluctuations of a+ for the mode k where their spectrum is 
maximum are plotted as the bold line at the Figure ^. Notice an enormous expo- 
nential growth of the fluctuations within a single background oscillation. Indeed, 
the amplitude "P+(fc) increases by factor 10 10 ! There are two factors which con- 
tributes to such a strong instability of fluctuations in the model. First, oscillating 
background fields are crossing the region with significant negative curvature of 
the effective potential, which results in tachyonic instability. Second, this region 
turn out to be a turning point for the background oscillations, where the fields 
spend significant portion of the oscillation. As a result, tachyonic instability is 
lasting long enough to make the backreaction of the fluctuations to be signif- 
icant already within single background oscillation. The regime of background 
oscillations will even not be settled. Therefore practically from the beginning we 
have to use the lattice simulations to study nonlinear dynamics of the fields. The 
results of full lattice simulations in the model are plotted in Figure |[ The simula- 
tions showed that the homogeneous fields 4> and er+ initially followed the classical 
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trajectory (Q) but, within one oscillation of the inflaton field, fluctuations grew 
too large to speak meaningfully of the fields as homogeneous oscillators. These 
fluctuations grew in such a way that c+ = <t_* almost exactly throughout the 
simulation. In other words Re 5a+ and Im (5<r_ were excited while Im <5cr+ and 
Re 5a- were not. Because of this we only plot the fields 4> an d cr+. 

In jjj we develop a general theory of tachyonic preheating, which occurs due 
to tachyonic instability in the theories with spontaneous symmetry breaking. 
Our approach combines analytical estimates with lattice simulations taking into 
account all backreaction effects. The process of spontaneous symmetry breaking 
involves transfer of the potential energy into the energy of fluctuations produced 
due to the tachyonic instability. We show that this process is extremely efficient 
and requires just a single oscillation of the scalar field falling from the top of the 
effective potential. We considered preheating in the hybrid inflation scenario, 
including SUSY motivated F-term and D-term inflationary models. 



3 Development of Equilibrium 

The character of preheating may vary from model to model, e.g. parametric ex- 
citation in chaotic inflation jjl and tachyonic preheating in hybrid inflation Q, 
but its distinct feature remains the same: rapid amplification of one or more 
bosonic fields to exponentially large occupation numbers. This amplification is 
eventually shut down by backreaction of the produced fluctuations. The end 
result of the process is a turbulent medium of coupled, inhomogeneous, classi- 
cal waves far from equilibrium. Despite the development of our understanding 
of preheating after inflation, the transition from this stage to a hot Friedmann 
universe in thermal equilibrium has remained relatively poorly understood. The 
details of this thermalization stage depend on the constituents of the fundamen- 
tal Lagrangian (P) and their couplings, so at first glance it would seem that a 
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Fig. 3. Means and variances in units of cf> c . The squared means {<f>} 2 and (cr+) 2 are 
ordinary solid lines while the field variances are thick lines. The mean of <j) starts at <f> c , 
oscillates once, and then decays. The mean of cr+ grows in antiphase to <f) and freeses 
at (j> c ■ 
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description of this process would have to be strongly model-dependent. Recently 
we performed a fully nonlinear study of the development of equilibrium after pre- 
heating . We have performed lattice simulations of the evolution of interacting 
scalar fields during and after preheating for a variety of inflationary models. We 
have found, however, that many features of this stage seem to hold generically 
across a wide spectrum of models. Indeed, at the end of preheating and beginning 
of the turbulent stage , the fields are out of equilibrium. We have examined 
many models and found that at there is not much trace of the linear stage of 
preheating and conditions at t* are not qualitatively sensitive to the details of 
inflation. We therefore expect that this second, highly nonlinear, turbulent stage 
of preheating may exhibit some universal, model-independent features. Although 
a realistic model would include one or more Higgs- Yang-Mills sectors, we treat 
the simpler case of interacting scalars. 

We have numerically investigated the processes of preheating and thermal- 
ization in a variety of models and determined a set of rules that seem to hold 
generically. These rules can be formulated as follows (in this section we use 
notations (f> = @i for the inflaton field and Xi °~ f° r other scalars <Pi) 

1. In many, if not all viable models of inflation there exists a mechanism for 
exponentially amplifying fluctuations of at least one field X- These mechanisms 
tend to excite long-wavelength excitations, giving rise to a highly infrared spec- 
trum. 

The mechanism of parametric resonance in single-field models of inflation 
has been studied for a number of years. This effect is quite robust. Adding 
additional fields (e.g. a fields) or self-couplings (e.g. x 4 ) nas little or no effect 
on the resonant period. Moreover, in many hybrid models a similar effect occurs 
due to tachyonic instability. The qualitative features of the fields arising from 
these processes seem to be largely independent of the details of inflation or the 
mechanisms used to produce the fields. 

2. Exciting one field x *s sufficient to rapidly drag all other light fields with which 
X interacts into a similarly excited state. 

We have seen this effect when multiple fields arc coupled directly to x an d 
when chains of fields are coupled indirectly to x- AH it takes is one field being 
excited to rapidly amplify an entire sector of interacting fields. These second 
generation amplified fields will inherit the basic features of the x field, i.e. they 
will have spectra with more energy in the infrared than would be expected for 
a thermal distribution. 

3. The excited fields will be grouped into subsets with identical characteristics 
(spectra, occupation numbers, effective temperatures) depending on the coupling 
strengths. 

We have seen this effect in a variety of models. For example in the models 
which we are going to consider the x an d °~ fields formed such a group. In 
general, fields that are interacting in a group such as this will thermalize much 
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more quickly than other fields, presumably because they have more potential to 
interact and scatter particles into high momentum states. 

4- Once the fields are amplified, they will approach thermal equilibrium by scat- 
tering energy into higher momentum modes. 

This process of thermalization involves a slow redistribution of the particle 
occupation number as low momentum particles are scattered and combined into 
higher momentum modes. The result of this scattering is to decrease the tilt 
of the infrared portion of the spectrum and increase the ultraviolet cutoff of 
the spectrum. Within each field group the evolution proceeds identically for all 
fields, but different groups can thermalize at very different rates. 

Here we will illustrate these results with a simple chaotic inflation model with 
a quartic inflaton potential. The inflaton <j) has a four-legs coupling to another 
scalar field \, which in turn can couple to two other scalars o\ and 02. The 
potential for this model is 

V = + i 9 Vx 2 + \h\x 2 °l + \h\x 2 o\ (9) 

Preheating in this theory in the absence of the Oi fields is well studied. For 

2 

nonsmall ^- the field % will experience parametric amplification, rapidly rising 
to exponentially large occupation numbers. In the absence of the \ ne ld (or 
for sufficiently small g) <j> will be resonantly amplified through its own self- 
interaction, but this self- amplification is much less efficient than the two-field 
interaction. The results shown here are for A = 9 x I0~ 14 (for CMB anisotropy 
normalization) and g 2 = 200A. When we add a third field we use h\ = lOQg 2 
and when we add a fourth field we use h\ = 200g 2 . 

One of the most interesting variable to calculate is the (comoving) number 
density of particles of the fields n(t) and their occupation number rife. The evo- 
lution of the total number density of all particles n to t is an indication of the 
physical processes taking place. In the weak interaction limit the scattering of 
classical waves via the interaction \g 2 $ 2 x 2 can be treated using a perturbation 
expansion with respect to g 2 . The leading four- legs diagrams for this interaction 
corresponds to a two-particle collision (cf)x — > cf>x), which conserves ntot- The 
regime where such interactions dominate corresponds to "weak turbulence" in 
the terminology of the theory of wave turbulence. If we see n to t conserved it 
will be an indication that these two-particle collisions constitute the dominant 
interaction. Conversely, violation of n to t(t) = const will indicate the presence of 
strong turbulence, i.e. the importance of many-particle collisions. Such higher 
order interactions may be significant despite the smallncss of the coupling pa- 
rameter g 2 (and others) because of the large occupation numbers rife. Later, 
when these occupation numbers are reduced by rescattering, the two-particle 
collision should become dominant and n to t should be conserved. For a bosonic 
field in thermal equilibrium with a temperature T and a chemical potential /z 
the spectrum of occupation numbers in the limit of classical waves is given by 
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Fig. 4. Time evolution of number density of particles in the model (|9|). The curves 
represent n,/,, n x , n CT1 , n CT2 from top to bottom. Unit of (conformal) time is a ■ 10 -36 
sec. 



Figure ^ shows an exponential increase of n(t) during preheating, followed 
by a gradual decrease that asymptotically slows down. This exponential increase 
is a consequence of explosive particle production due to parametric resonance. 
After preheating the fields enter a turbulent regime. In our simulations we see 
n(t) decreasing during this stage. This decrease is a consequence of the many- 
particle interactions beyond the four legs rescattering. 

Figure || illustrates the late-time spectrum of particles. Initially spectra are 
concentrated at the infrared section, which is gradually flattening as it ap- 
proaches a thermal distribution ( |ic| ) (i.e. a slope of —1 to depending on the 
chemical potential and the mass). The spectra of three fields are essentially 
identical, which leads to 

n x w n al w n a2 < n^. (11) 



Occupation Number: t-2400. 
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Fig. 5. Occupation number versus momentum for the model ([]). The curves represent 
spectra of particles for the fields <f>, \, Ci and 02 from top to bottom All the fields other 
than the inflaton have nearly identical spectra 
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Another important point is that the interaction of \ an d &i does not affect the 
preheating of xu but does drag <ii exponentially quickly into an excited state. 
The fields 0$ are exponentially amplified not by parametric resonance, but by 
their stimulated interactions with the amplified x field. Unlike amplification by 
preheating, this direct decay nearly conserves particle number, with the result 
that n x decreases as Oi grow. 

Interacting waves of scalar fields constitute a dynamical system. Dynamical 
chaos is one of the features of wave turbulence. In || we address the question 
how and when the onset of chaos takes place after preheating. To investigate the 
onset of chaos we have to follow the time evolution of two initially nearby points 
in the phase space. Consider the theory with the potential (^|) with two fields 
4> and x only (which we collectively denote as /). Consider two configurations 
of a scalar field / and /' that are identical except for a small difference of the 
fields at a set of points xa- Chaos can be defined as the tendency of such nearby 
configurations in phase space to diverge exponentially over time. This divergence 
is parametrized by the Lyapunov exponent for the system, defined as 



A = -log — — 
t y A 



(12) 



where A is a distance between two configurations and Aq is the initial distance 
at time 0. In the context of preheating is is convenient to define A as 



f_A ~ fA 
f' A + fA 



I'a ~ h 
f'A + h. 



(13) 



that is well regularized even while the field x is being amplified exponentially 
during preheating. Figure || shows the Lyapunov exponent A. We see the onset 
of chaos only at the end of preheating. The plot for the 4> field is nearly identical. 



The Lyapunov exponents for the fields were A^ 
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Fig. 6. The Lyapunov exponent A for the fields cj> and x using the normalized distance 
function A. 
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adopted in the simulation). This corresponds to a very fast onset of chaos around 
the moment t* where the strong tubulence begins. 

The highlights of our study for early universe phenomenology are the fol- 
lowing. The mechanism of preheating after inflation is rather robust and works 
for many different systems of interacting scalars. There is a stage of turbulent 
classical waves where the initial conditions for preheating are erased. Initially 
before all the fields have settled into equilibrium with a uniform temperature, 
the reheating temperature may be different in different subgroups of fields. The 
nature of these groupings is determined by the coupling strengths. 

4 Parametric Excitations of Fermions 

The interaction between bosons and inflaton(s) may results in the copious pro- 
duction - preheating - of the bose particles, either due to the parametric res- 
onance as in the chaotic inflation or due to the tachyonic instability as in the 
hybrid inflation. Preheating of bosons has been studied in details with analytic 
approximations as well as with the lattice simulations. Consider second line of 
Eq. (]]]) and ignore mixing between fermions x an( i gravitino ip^ (which cor- 
responds to the rigid SUSY limit). There is interaction between inflaton and 
fermions. Let us consider a simple Yukawa interaction hiptptp, which should lead 
to the production of fermions from inflaton oscillations (in this section we use 
notation <fi for the inflaton field). For fermions, the Pauli exclusion principle 
prohibits the occupation number from exceeding 1. For this reason, it has been 
silently assumed that fermions are created in the three-legs perturbative process 
<j) — > XX where individual inflatons decay independently into pairs of x-particlcs. 
Following (^] let us, however, consider the Dirac equation for a massless quantum 
Fermi field x(t, x) infracting with the background inflaton 



where V M is the derivative with the spin connection. Here, similar to the bosonic 
case, the inflatons producing fermions also act not as individual particles but as a 
coherently oscillating field 4>(t). Indeed, the equation for the eigenfunctions of the 
quantum fluctuations in this theory can be reduced to a second-order equation 
for an auxiliary field X(t, x), so that \P = [7 M V M + h<p] X. The eigenmodes of 
the auxiliary field have the form Xk(t)e +lk x i? r , where the R r are eigenvectors 
of the Dirac matrix 7 with eigenvalue +1 (we are using the representation 
of gamma matrices where 70 = diag(i, i, — i, — i)). The temporal part of the 
eigenmode obeys an oscillator-like equation with a complex frequency which 
depends periodically on time 



[7"V M + H{t)\ x = , 



(14) 




(15) 



where fl = y/qf. The background oscillations enter in the form f(t) — 
where 4>o is their amplitude. Without imaginary part, this equation will the the 
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Fig. 7. The occupation number nk of fermions in rn^cj) 2 -inflation. The initial resonance 
parameter is qo = 10 6 and the mode is kq ~ y^ 771-2 ' ^ = "m" ■ 

equation for the fluctuations of the bose field, which results in its parametric 
resonance. The imaginary part of the frequency in Eq. (|l^) guarantees the Pauli 
blocking for the occupation number nu ■ The results for rik can be formulated as 
follows H . Even though the Yukawa interaction contains a small factor h, one 
cannot use the perturbation expansion in h. The strength of the effect is ulti- 
mately determined by the dimensionless parameter q. For instance, for chaotic 
inflation with A</> 4 potential, q — h- is typically very large. For instance, in the 
simple supergravity theory with W = V~\<j) 3 and single chiral multiplet, h = ^/2\ 
and q = 2. For the quadratic inflaton potential m 2 (j) 2 , the parameter q = ^j 2 -, 
where </>o is the amplitude of the inflaton oscillations, which initially is ~ O.lMp. 
The modes get fully excited with occupation numbers ~ 1 within tens of 
oscillations of the field 4>, and the width of the parametric excitation of fermions 
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K 2 

Fig. 8. The comoving occupation number of fermions in m^cf) 2 -inflation after 50 in- 
flaton oscillations for initial resonance parameter q = 10 3 . Expansion of the universe 
destroys the details of the resonance band and leads to a fermi- sphere of width q 1 ^ 4 m. 
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in momentum space is about g 1 / 4 m e /j, where m e ff is the frequency of the back- 
ground inflaton oscillations. It turns out that the analytic theory of fermionic 
preheating can be advanced similar to that for the bosonic preheating, see || 
for details. For illustration in Figure we show the time evolution of the occu- 
pation number of fermions in m 2 <fi 2 inflation, where excitations of fermions is a 
stochastic process. Figure ^ shows the spectrum of fermionic excitations, which 
stochastically filling a (Fermi) sphere in the momentum space . 



5 Gravitino from Preheating 

Let us now consider the third line in the supergravity Lagrangian d), which 
describes the gravitino field ip^. Assume that the gravitino mass m 3 / 2 ~ 10 2 — 
10 3 Gev. Such particles decay very late after nucleosynthesis and lead to disas- 
trous cosmological consequences unless the abundance of gravitinos (in units of 
the number density of photons) is extremely suppressed, n^^/n^ < 10 -15 . This 
puts constrains on the reheat temperature Tr to avoid gravitino overproduction 
from thermal scattering of particles in the very early universe, Tr < 10 8 Gev. 
However, gravitinos can be produced not only in the thermal bath after preheat- 
ing but even earlier during the inflaton oscillations, similarly to the production of 
fermions from preheating |||4| . Thus the investigation of the non-thermal grav- 
itino production in the early universe may serve as a useful tool to discriminate 
among various versions of inflations. 

To study gravitino production from preheating, first one shall derive its equa- 
tion, which is a challenge by itself. The formalism for gravitinos in an expanding 
universe and in the presence of complex scalar fields with non-vanishing VEV's 
was recently formulated in paper || . The equation for the gravitino has on the 
left hand side the kinetic part = e^ v P<T 75 7„ £> p Vv, and a rather lengthy right 
hand side (we use the long derivative 2? M with the spin connection and Christof- 
fel symbols, for which T>^ v — 0). Apart of varying gravitino mass m.3/2 = ^gr, 
the right hand side contains a chiral connection and various mixing terms like 
those in the 3rd line of (|l|). For a self-consistent setting of the problem, the 
gravitino equation should be supplemented by the equations for the fermions 
Xi and gauginos A Q mixing with gravitino, as well as by the equations deter- 
mining the gravitational background and the evolution of the scalar fields. We 
formulated cosmological extension of the super-Higgs effect Q . It turns out that 
supersymmetry is spontaneously broken if the combination 

H 2 + m 2 3/2 > (16) 

is positive. In flat space-time usually supersymmetry breaking is associated with 
the non- vanishing gravitino mass TO3/2. In an expanding universe the Hubble pa- 
rameter H plays an equally important role. For instance, the gravitino mass (i.e. 
superpotential W) can vanish, but supersymmetry will be still spontaneously 
broken. For further discussion of the nonthermal gravitino production, let us 
make some simplifications. We consider the supergravity multiplet and a single 
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chiral multiplet containing a complex scalar field 4> with a single chiral fermion 
X- This is a simple non-trivial extension which allows us to study the gravitino 
with a non-trivial FRW cosmological metric supported by the scalar field. A 
nice feature of this model is that the chiral fermion \ can be gauged to zero 
so that the mixing between ip^ and \ m (0) is absent. We also can choose the 
non-vanishing VEV of the scalar field (inflaton) in the real direction. First we 
will derive the equation for a spin 3/2 field in a curved background metric with 
non-vanishing VEVs for the scalar fields. From (jl]) we can obtain the equation 
for the gravitino We use a plane- wave ansatz ?/> M ~ e lk x for the space-dependent 
part. Component 4>o is related to ipi via the constraint equation. Next, ipi can be 
decomposed into its transverse part ipf, and to the longitudinal part x/jf which 
is defined by the trace 7 l, 0i- Two degrees of freedom of ijj^ are associated with 
the transverse part ipf , which correspond to helicity ±3/2 and two degrees of 
freedom are associated with Y^Pi ( or "00 ) which correspond to helicity ±1/2. 
Equation for the helicity ±3/2 states for /i = i is 

(l^ + ^° + m 3/2 a^f = 0. (17) 

The transformation ijjf = a -1 / 2 !?^ reduces the equation for the transverse part 
to the free Dirac equation with a time-varying mass term m 3 / 2 a, c.f. eq. (|l4|). 
In the previous section we explained how to treat this type of equation. The 
essential part of is given by the time-dependent part of the eigenmode of the 
transversal component _Xy(t), which obeys second-order equation (c.f. (JlF 



X T + (fc 2 + (m 3/2 a) 2 - i(m 3 / a a)) X T = . (18) 

Production of helicity 3/2 gravitino is essential at the fast rolling stage of inflaton 
evolution, when W may be changed nonadiabatically. 

The corresponding equation for gravitino with helicity 1/2 is more compli- 
cated. It is convenient to use combination 9 = 7^ which fully describes helicity 
1/2 states. Equation for combination The equation for 9 is 



d t + S-ik-7 7o A) 6» = 0, (19) 

where B = -&A- ^70 (1 + 3 A), and 

, p-3m 2 Ml 2mM 2 p 

A = 2m + 7o ^ q £j2 = A i + 70^2 ■ 20 

p + 3m^Mp p + 6m z Mp 

Here p and p are the energy-density and pressure of the background scalar field. 
For the single scalar field \A\ 2 = 1 for an arbitrary superpotential W . Thus 

A can be represented as A — — exp ^2i dt fJ>(rf)j , where /i — WW + A, 

the correction A — 0(Mp 1 ) is given in g. The time-dependent factor of the 
spinor 9, which we denote as fk(t), obeys a second-order differential equation. By 



Preheating, Thermalization and Supergravity 15 



the substitution f k (t) = E(t)X(t), with E = {-A*) 1 / 2 e I dtReB ^ the equation 
for the function fk(t) is reduced to the oscillator-like equation for the time- 
dependent mode function Xl of the longitudinal component: 

x L + (k 2 + n\- ifi L ) x L = o (2i) 



= /i— |i? sin 2 J fidt—im (l + 3 cos 2 j fidi). Equation ( |2l| ) describes 
the creation of gravitinos from preheating. Notable the gravitino with helicity 
1/2 remains coupled to the changing background even in the limit Mp — > oo. 
In a sense, the gravitino with helicity 1/2 remembers its goldstino nature. Thus 
gravitino production in this background in general is not suppressed by the 
gravitational coupling. The main dynamical quantity which is responsible for 
the gravitino production in this scenario will not be the small changing gravitino 
mass TO3/2, but the mass of the chiral multiplet fi, which is much larger than 

But do we really produce gravitinos from preheating, or we describe the 
production of the chiral fermion? To answer this question, we have to investigate 
more realistic and more complicated problem with several chiral multiplets. As 
the first step in this direction we derived equations for gravitinos for an arbitrary 
number of chiral multiplets plus SYM sector j^j . For instance, generalization of 



the equation (19) for the next simplest case of two chiral multiplets even in the 



limit Mp — > 00 is reduced to the complicated second-order spinor equation 

%+B-tis.- 77 i) 9 + G \dt + (k 2 - ik ■ 77 i)j 9 = 0, (22) 



the matrix G is constructed from the backgroud scalars. We expect solutions 
of this equations will allow us to find the residual gravitino production. Recent 
report |g| indicates that helicity 1/2 gravitinos are produced in the model with 
two chiral multiplets. 
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